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Abstract 

The recently developed chiral nucleon-nucleon (NN) potential by Epelbaum et al. has been 
employed to study the two-nucleon bound and scattering states. Chiral NN potential up to next- 
to-next-to-next-to leading order (N 3 LO) is used to calculated the np differential cross section and 
deuteron binding energy in a realistic three dimensional approach. The obtained results based 
on this helicity representation are compared to the standard partial wave (PW) results. This 
comparison shows that the 3D approach provides the same accuracy in the description of NN 
observables and the results are in close agreement with available experimental data. 
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I. INTRODUCTION 



It is a long time that the standard PW decomposition has been used to solve the few- 
body problems. In this approach one should sum all PW's to infinite order, but in practice 
one truncates the sum to a finite angular momentum number which is dependent to the 
considered energy. It means that in higher energies one will need many PW components, 
which leads to very complicated expressions, to achieve the convergence results. It appears 
therefore natural to avoid the very involved angular momentum algebra which is inherent in 
the PW representation of permutations, transformations and especially the 3N forces and 
work directly with vector variables To this aim in the past decade the main steps have 
been taken by Ohio-Bochum collaboration (Elster, Glockle et al.) and Bayegan et al. to 
implement the 3D approach in few-body bound and scattering calculations (see for examples 
Refs. {^-{sjl). The 3D approach replaces the discrete angular momentum quantum numbers 
with continuous angle variables and consequently it takes into account automatically all 
PWs. So in contrast to the truncated PW approach, the number of equations in the non 
truncated 3D representation is energy independent. Therefore this non PW method is more 
efficient and applicable to the three- and four-nucleon scattering problems which consider 
higher energies than the corresponding bound state problems. It should be clear that the 
building blocks to the few-body calculations without angular momentum decomposition are 
two-body off-shell t-matrices, which depend on the magnitudes of the initial and final Jacobi 
momenta and the angle between them. Fachruddin et al. have formulated the NN bound 
and scattering states in a 3D representation and they have numerically illustrated the np 
differential cross section and deuteron binding energy by using two realistic model interac- 
tions, i.e. the Bonn-B and the AV18 

aa 

. They have incorporated the momentum vectors 
directly into the bound and scattering equations, and the total spin of the two nucleons 
is treated in a helicity representation with respect to the relative momentum of the two 
nucleons. Despite many successes that conventional approaches achieved in incorporating 
the NN potentials, like the CD-Bonn, the Nijmegen I and II, and the AV18, in nuclear 
structure and reaction calculations, there are certain deficiencies that require a reliable ap- 
proach which is based on the theory of strong interactions, the quantum chromodynamics 
(QCD). These deficiencies can be categorized as no connection to QCD, model-dependent 
with the lack of 3iV force to be add on, the gauge and chiral symmetries hard to be reached 
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and finally fine tuning in not achievable order by order of increasing momenta. Based on 
the spontaneously and explicitly broken chiral symmetry it is possible to construct nuclear 
forces in the 
by Weinberg 



et al 



■ H, 



ramework of the chiral perturbation theory. This approach has been founded 



10. 



ll| and further expanded 



jy 



and recently by Epelbaum et al. 
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Ordonez et al. jl2|], Kaiser et al. [l3|, Entem 



161 ]. In order to compare the 3D and the PW 



approaches in a more fundamental basis, we intend to incorporate the new chiral potential 
|l7j ] into the 3D few-body calculations. In the first step we are preparing this potential in 
an appropriate operator form, which is consistent with 3D representation, to calculate the 
np differential cross section and also the deuteron binding energy. 

II. A BRIEF REVIEW OF THE 3D FORMALISM FOR NN BOUND AND SCAT- 
TERING STATES 



The NN differential cross section is given as: 



^ = ( 27r )4(™)2 
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a (p' m Sl m S2 m' ti m' t2 \T\pm Sl m S2 m tl m t2 ) a 



(1) 



where m Si and m ti indicate the projection of the spin and isospin of the nucleons, p and 
p' are initial and final relative momentum of the two nucleons and the operator T is the 
2N transition matrix determined by Lippmann-Schwinger equation. In order to calculate 
the NN differential cross section we need to calculate the matrix elements of the physical 
representation of NN T-matrix, i.e.: 



a (p' m Sl m S2 m' ti m' t2 \T\pm Sl m S2 m tl m t2 ) a , 



which are given in the antisymmetrized basis states, i.e. |p m Sl m S2 m tl mt 2 ) t 



(2) 
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Pi2)\p m Sl m S2 m tl m t2 ) . These matrix elements can be obtained by a summation over the 
on-shell momentum helicity T-matrices multiplied with the rotational matrices and Clebsch- 

n / n 

Gordon coefficients ]4]. As indicated in Ref. |3J the projection of the Schrodinger equation 
on the helicity basis states leads to the coupled integral equations which, after simplification, 
are actually only two dimensional integral equations. Since the deuteron has spin 1 there 
are three possible values for the helicity projections, namely A = — 1, 0, +1. The symmetry 



properties allow to consider only A = +1, 0. Thus one obtains a set of two coupled integral 
equations in two variables, the magnitude of the relative momentum vector, i.e. p, and the 
angle between p and the arbitrarily chosen z-axis, i.e. 9. 



III. PREPARATION OF THE CHIRAL POTENTIAL IN 3D REPRESENTATION 

The general form of the NN potential by considering the rotation, parity, and time 
reversal invariance can be written as a linear combination of six f2j operators, which are 
consistent with the helicity basis representation j^]: 

6 

(p'|V|p> = V(p',p) = $>(p',p,7) X) (3) 

i=i j 

where Vi(p',p, 7) are scalar functions which depend on the magnitudes of p and p' and 
also the angle between them, 7 = p.p' , and A is a 6 x 6 matrix. The Qi operators are: 

^ = 1, q 2 = s 2 , q 3 = s • p's ■ p', n A = s • p's ■ p, n 5 = (s ■ p') 2 (s ■ p) 2 , q 6 = s ■ ps ■ p. By 

this representation the spin-dependent parts of the matrix elements of the potential can be 
easily evaluated in the helicity basis states. We intend to use the chiral AW potential up 
to N 3 LO of chiral expansion which consists of "one- and two-pion exchanges (1PE, 2PE) 



and a string of contact interactions with an increasing number of derivatives 



zero, two, 



four) that parameterize the shorter ranged components of the nuclear force" [17] . In order 
to use the chiral potential in 3D formalism we need first to rewrite this potential in an 
appropriate operator form which is consistent with helicity representation. To this aim we 
should overcome the following two possible issues: 

• the calculation of original low energy coefficients (LEC's) for incorporating the contact 
terms 

• the representation of the spin dependent parts in term of fij operators 
As indicated in reference 



171 | the chiral potential at N 3 LO consists of contact terms 
which contain 24 original LEC's; C5, CV, Ci, Cj and D±, D 15 . In order to calculate 
any observable with chiral potential in the PW approach it is sufficient to project only the 
contact interactions in the fourteen PW channels up to J = 3 and it is not necessary to 
consider the higher channels, i.e. J = 4,5,6, etc. Once the 24 spectroscopic LEC's have 
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been determined by fitting to the phase shifts of the Nijmegen potential, the original ones 
can be obtained uniquely. This is a serious problem to apply the contact terms of the 
chiral potential in 3D approach, since in this approach we consider all of the PW channels 
simultaneously. To overcome this problem we have used the connection between the 3D and 
PW representations of matrix elements of the iViV potential [18J to sum over these fourteen 
channels and to obtain the matrix elements of the potential in momentum helicity basis. 

In order to make the chiral potential compatible with helicity representation, we highlight 
the spin dependent parts of the potential as follows: 

cr\ ■ ct>- cr 1 ■ q cr 2 ■ q, cr 1 ■ k <x 2 • k, 

i((T 1 + cr 2 ) ■ (q x k), ct ± ■ (q x k) cr 2 ■ (q x k). (4) 

These parts can be simply written in term of Qi operators as: 

(Ti ■ cr 2 = 2Q 2 — 3fii, (5) 

p'pa 2 ^ 2gV7-p) 
cri • q<x 2 • q = -blli H "2 H "3 

7 7 

-2ppll 4 -\ il 5 H IIq, (6) 

7 7 

c p'pa 2 p'{p'i + p) 
at ■ kcr 2 • k = -7O1 fi 2 H 3 

4 47 27 

+^^ 4 - P ^n 5 + P^ + ^n,, (7) 

2 2 7 5 2 7 6 ' v ; 

i(<Ti + cr 2 ) • (q x k) = z(<7! + <r 2 ) • (p x p) 

p'pa 2 2p'p , 2p'p 2p'p 
= ih il 3 + 2p pill -\ \l 5 H il 6 , (8) 

7 7 7 7 

<Ti ■ (q x k) cr 2 ■ (q x k) = cri ■ (p x p) cr 2 • (p x p) 

= -p'Va^i + P 2 p 2 a 2 n 2 + 2p' 2 p 2 jQ 4 - 2p /2 p 2 n 5 , (9) 

where 7 = p' • p, a = a/1 — 7 2 , 6 = p /2 + — 2p'p7 = q 2 , c = p' 2 + p 2 + 2p'p , y = 4k 2 . 
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IV. NUMERICAL RESULTS 



In this section we present the obtained numerical results for np differential cross section 
and deuteron binding energy with chiral potential up to N 3 LO in the 3D approach. In order 
to demonstrate the effectiveness of 3D formalism we have compared our numerical results 
with the corresponding PW results as well as the experimental data. The low energy co- 
efficients in the chiral potential are determined for given cut-off parameters Aj and A 2 by 
fitting to NN data, where the cut-off Ai regulates the high-momentum components of the 
interacting nucleons and the cut-off A 2 which appears in the spectral function regulariza- 
tion excludes the high-momentum components of the two-pion exchange. We use in our 
calculations different combinations of (A 1; A 2 ) for N 3 LO as: (450,500), (600,600), (550,600), 
(450,700), (600,700) in units of MeV/c. 

TABLE I: Deuteron binding energy calculated for the chiral potential at N 3 LO in 3D approach for 
three different cut-off sets in comparison with the PW and experimental results. 



(Ai,A 2 ) 


3D [MeV] 




PW [19] [MeV] 


(450,500) 


-2.216 




-2.215 


(450,700) 


-2.219 




-2.218 


(600,700) 


-2.222 




-2.220 


PW [17] 
EXP 




(-2.216X-2-223) 
-2.224575(9) 





The numerical results for np differential cross section in four different energies of the 
projectile in the laboratory system and for different cut-off sets have been shown in figure 
([!]). In the first row of this figure we have presented a comparison between 3D and PW 
results for Ei a t, = 50 MeV. The PW results have been taken from [l^] where the calculation 
is up to Jmax = 6. Both 3D and PW results are in good agreement in forward and backward 
angles. Also the comparison of both approaches with experimental data in backward angles 
shows a very close agreement. In the next rows the same comparison has been shown for 
Eiab = 96, 143 and 200 MeV, and as we see both approaches match together and also to 
the experimental data perfectly. In table [J we have presented our numerical results for 
deuteron binding energy in comparison with PW and experimental data. Our numerical 
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results for three cut-off sets with the values —2.216, —2.219 and —2.222 MeV are in good 
agreement with the very recent corresponding PW results [3], and also with Epelbaum et 



al. PW achievements 



171 ] . The agreement between the 3D and PW results as well as the 



experimental data is quiet satisfactory. 

Although we have studied the 2N systems we conclude that the 3D approach is promising 
to be simpler for more complex few-body systems by providing a strictly finite number of 
coupled three-dimensional integral equations to be solved. The number of the integral 
equations in the 3D approach is consistent and do not depend on the energy of the system. 
This subject is more important when we consider 3N and AN scattering problems in which 
the number of the equations in the higher energies makes the problems more complex. The 
recently developed 3iV bound state in 3D approach [t| can be used to calculate the 3 H and 
3 He binding energies by using the chiral potential. The 3N scattering, Nd capture and 3N 
photodisintegration calculations with this new form of the chiral potential are interesting 
goals that we are pursuing. Also the incorporation of 3N chiral forces in novel 3D approach 
calculations is one of other interesting problems that can be done. 
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FIG. 1: The np differential cross section for four different energies. The left figures are 3D and 
the right ones are PW results. In each figure the 3D or PW numerical results are obtained for 



five cut-off sets, i.e. Cl:(450,500), C2:(600,600), C3:(550,600), C4:(450,700), C5:(600,700) 
The PW results have been obtained up to J = 6 and they have been taken from ref. 
experimental data have been taken from ref. [2o| for Ei a b = 50 MeV (EXP1), from refs 
for Ei ab = 96 MeV (EXP 1 and EXP 2 respectively), from ref. 



and from refs. 



24, 



21 



MeV. 
he 



22] 



23(] for E iab = 142.8 MeV (EXP1) 



25[ for Eiab = 200 MeV (EXP1 and EXP 2 respectively). 
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